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ABSTRACT 

We present theoretical results for the expected fraction of Lyman Break Galaxies 
(LBGs) to be detected as strong Lyman-a emitters (LAEs) in the redshift range 
5 < z < 7. We base our analysis on the 2-billion particle SPH simulation MareNostrum 
High-z Universe. We approximate galaxies as static dusty slabs with an additional 
clumpy dust distribution affecting stellar populations younger than 25 Myr. The model 
for the Lyman-a escape fraction is based on the results of our Monte-Carlo radiative 
transfer code (CLARA) for a slab configuration. We also fix the transmission of Lyman-a 
photons through the intergalactic medium to a constant value of 50% at all redshifts. 
From the results of this model we calculate ccLya , the fraction of Lyman Break Galaxies 
with Lya equivalent width (EW) larger than 50A. We find a remarkable agreement 
with observational data at 4.5 < z < 6. For bright (—22 < Mtjy < —20.5) and 
faint (—20.5 < Muv < —18.5) galaxies our model predicts x^ya = 0.02 ± 0.01 and 
Xhya = 0.47 ± 0.01 while observers report x^ya = 0.08 ± 0.02 and XLya — 0.47 ± 0.16, 
respectively. Additional evolution of the extinction model at redshift z ^ 7, that 
decreases the intensity of transmitted Lyman-a radiation by a factor of /t = 0.4 as to 
match the LAE luminosity function at z ~ 6.5, naturally provides a good match for 
the recently reported XLya fractions at z > 6.3. Exploring different toy models for the 
Lyman-a escape fraction, we show that a decreasing Lyman-a escape fraction with 
increasing UV galaxy luminosity is a key element in our model to explain the trend 
of larger XLya fractions for fainter LBGs. 

Key words: galaxies: high-redshift — galaxies: evolution — methods: N-body sim- 
ulations 



1 INTRODUCTION 

The study of distant star forming galaxies is being driven by 
even more sensitive observations. Observational samples are 
now commonly gathered both for Lyman Break Galaxies 
(LBG) and Lyman-a e mitting (LAE) gala xies in the red- 
shift range z > 4 (see iHaves et al.l |20lJ) and references 
therein). Nevertheless, a clear physical connection between 
these two populations is difficult to establish given the com- 
plex physics involved in the trans mission Lyman-a rad iation 
in the interstellar medium (ISM) (jHansen fc Qhll2006h . 

Recently, observational results have reported on the 
fraction of LBGs that show strong Lyman-a emission, xi^ya, 
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at a given absolut e rest-frame U V magnitude (jStark et al.l 
I2010l .i2011: Schcn kor et al.ll201lf ). This fraction derived from 
observational data provides a simple test to theoretical mod- 
els that seek to explain, under a sound physical model, the 
connection between LBGs and LAEs. Because this test is 
concerned with fractions of a population, it is largely insensi- 
tive to changes in absolute number densities of the observed 
galaxies, adding new information with respect to analysis 
based on the LAE luminosity function (LF). 

In this paper we make theoretical predictions on the 
XLya fraction as function of absolute rest-frame magni- 
tude Muv in the redshift range 5 < z < 7. These re- 
sults are based on a simulation of a cosmological volume 
that follows the dynamical evolution of roughly 2 billion 
dark matter and gas particles including star formation and 
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supernovae feedback. The clumpy dust attenuation model 
and its effects on t he UV luminosity function were de- 
scribed in detail in iForero-Romero et alj (|201Ct ) (Paper I 
hereafter). The corresponding Lyman-a escape fraction is 
ca lculated using the results of our Monte-Carlo code CLARA. 
In iForero- Romero et al.l ()201lh (Paper II hereafter) we dis- 
cussed at length the structure of the code and the impli- 
cations for the luminosity functions in the redshift range 
5 < z < 7. 

In Section [2] we describe the numerical simulation and 
the model for the UV/Lyman-Q emission and its associated 
extinction/escape fraction. Next (^we describe the results 
of this model for the SLya statistic in the redshift range 
5 < z < 7. We also explore different toy models for the 
Lyman a escape fraction in order to better understand the 
trends found in the simulation. We discuss the results in 
Section |4] and present our conclusions in Section (5] 



2 A SIMULATION OF HIGH REDSHIFT LEGS 
AND LAES 

The cosmological simulation, the algorithm of galaxy finding 
and the spectral modeling (UV continuum and Lyman-a 
line) have been thoroughly described in Paper I and II. Here 
we summarize the most relevant features for this Paper. 



2.1 The MareNostrum High-Z Universe 
Simulation 

The MareNostrum High-z Universe simulatioi:Q follows the 
non linear evolution of structures in baryons (gas and stars) 
and dark matter, starting from z = 60 within a cube of 
50/i"^Mpc comoving on a side. The c osmological parame - 
ters used correspond to WMAPl data (|Spergel et al.ll2003l ') 
and are = 0.3, fib = 0.045, Qa = 0.7, erg = 0.9, a Hubble 
parameter h = 0.7, and a spectral index n = 1. The initial 
density field has been sampled by 1024'^ dark matter parti- 
cles with a mass of muM = 8.2 x lO'^/i"^M0and 1024^ SPH 
gas particles with a mass of mgas = 1.4 x lO^/i~^M0. The 
simulation has been perform ed using the TREEPM-I-SPH 
code GADGET-2 (|Springell2005l '). The gravitational smoothing 
scal e was set to 2 fe~^kpc in com oving coordinates. We fol- 
low ISpringel fc Hernguistl (|2003l ) to model the cooling, star 
formation and strong kinetic feedback model in the form of 
galactic winds. We identify the objec ts in the simulations 
using the AMIGA Halo Finder (AHF) (|Knollmann fc Kneb3 
I2OO9I ). All objects with more than 1000 particles, dark mat- 
ter, gas and stars combined, are used in our present analyses. 
We assume a galaxy is resolved if the object contains 200 or 
more star particles, which corresponds to objects with > 400 
particles of gas. 



2.2 UV and Lyman-a emission 

The photometric properties of galaxies are calculated em- 
ploying the stellar pop ulation synthesis model STARDUST 
( Devriendt et al.l 1993 ). Using the methods described in 
iHatton et all (|2003l '). 



We consider only the intrinsic Lyman-a emission asso- 
ciated to star formation. We assume that the number of Hy- 
drogen ionizing photons per unit time is 1.8 x 10^^ photons 
s~^ f or a star formation rate of 1 M0/yr (|Leitherer et al.l 
Il999h . 

Assuming that 2/3 of these photon s are converted t o 
Lyman-a photons (case-B recombination. [Osterbroc^l 19891 ) . 
the intrinsic Lyman-a luminosity as a function of the star 
formation rate is 



1.9 x 10*^ X (SFR/Mq yr~^) erg s"\ 



(1) 



A change in the escape fraction of ionizing photons can thus 
induce changes in the intrinsic Lyman-a emission. In our 
model, we take this ionizing photon escape fraction to be 
negligible. 



2.3 Dust attenuation 

Our approach to calculate the dust extinction is purely 
phenomenological. The extinction curve for each galaxy 
is different depending on its metallicit y and gas contents 
(|Guiderdoni fc Rocca-Volmerangdl 19871 ). The dust attenua- 
tion model parametrizes both the extinction in a homoge- 
neous ISM and in the molecular clouds around young stars, 
following the physical model of I Chariot fc Fall (|200cl ). The 
attenuation from dust in the homogeneous ISM assumes a 
slab geometry, while the additional attenuation for young 
stars is modeled using spherical symmetry. 

We first describe the optical depth for the homogeneous 
interstellar medium, denoted by t^^'^\\). We take the mean 
optical depth of a galactic disc at wavelength A to be 



(Nh) 



2.1 X 10^^ atoms cm 



where Ax/Av is the extinction curve from iMathis et al.l 
(|l983h . Zg is the gas metallicity, (Nh) is the mean atomic 
hydrogen column density and rj = {1 + z)^'^ is a, factor that 
takes into account the evolution of the dust to gas ratio at 
different redshifts. 

The extinction curve depends on the gas metallicity Zg 
and is based on an interpolation between the solar neighbor- 
hood and the Large and Small Magellanic Clouds (r = 1.35 
for A < 2OOOA and r = 1.6 for A > 2000A). 

Stars younger than a given age, tc, are subject to an 
additional attenuation in the birth clouds (BC) with optical 
depth 



rr(A)=(^-l)rr^^(A), 



(3) 
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where fj. is the fraction of the total optical depth for these 
young stars with respect to that is found in the homogeneous 
ISM. 



2.4 Lyman-a escape fraction 

In Paper II using our radiative transfer code CLARA we ob- 
tained the Lyman-a escape fraction for the correspond- 
ing slab described as a function of the product {aro)^^^^Ta, 
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where tq is the Hydrogen optical depth, Ta is the opti- 
cal depth of absorbing material (for albedo values of A, 
Ta — (1 — A)rd, where Td is the dust optical depth), and 
a is a measure of the gas temperature defined as a = 
Ai'L/(2Ai^_D), Avd ~ {vp/c)v(i is the Doppler frequency 
width, and Vp = {2kT /niHY^^ is \/2 times the velocity dis- 
persion of the Hydrogen atom, T is the gas temperature, 
vriH is the Hydrogen atom mass and Avl is the natural line 
width. 

We find that the equation 



where 



1 - exp(-P) 



1/3^ n3/4 



(4) 
(5) 



with e = 3.5 provides a reasonable description of the Monte- 
Carlo results for [aroY^^Ta < 200 for the slab geometry with 
homogeneously distributed sources. 

In addition to the foreground, homogeneous ISM ex- 
tinction, we model the attenuation due to the birth clouds. 
Under the physical conditions for these poor neutral hydro- 
gen clouds, T(f'^ ~ together with ar^^ < 1 and fur- 
thermore Ta'~' > 1, the enhancement in absorption by reso- 
nant scattering becomes irrelevant. Under these conditions, 
we take the escape fraction as the continuum extinction at 
A = 1260 A for a spherical geometry. 

The effect of gas kinematics is not taken into account 
in the model. Appropriate outfiow co nfigurations can in- 
crease the Lyman-g escape fraction (|Kunth et al.l 1 19981 : 
IVerhamme et al.ll2006l : lAtek et al.|[2008l ). 



2.5 Transmission through the IGM 

In Paper II we assumed that the IGM allows the transmis- 
sion of 50 per cent of the Lyman-a line at every redshift 
5 < z < 7. This is quantified by a transmission coefficient 
TiGM ~ 0.50, which can be considered optimistic given re- 
cent numerical estimations. 

Simplified analytic modeling jStark et al.l l201ll ) of 
TiGM ~ exp (tc) where T q is the optical depth for Lyman-a 
photons (calculated from lMeiksinI (|2006l )). yield Tigm = 0.58 
and Tigm = 0.51 at redshifts z ~ 5 and z ~ 6. 3D Monte 
Carlo radiative transfer calculations of individual galax- 
ies obtain values fo r the IGM transmis sion at 2 ~ 6 of 
Tigm = 0.26lo i8 (|Laursen et all l201ll '). a result that is 
readily explained by the effect of circumgalactic gas infa ll 
jPiikstra et all l2007l : llliev et all I2OO8I: iDaval et al.l I2OI j ). 
Radiative transfer calculations by IZheng et al. (|2010l ) find 
both a larger scatter and a lower mean value in the fraction 
of intrinsic Lyman-a radiation effectively detected at z ~ 6. 
However, in that calculation the authors need to boost by 
a factor of ~ 5 the intrinsic Lyman-a intensity in order 
to match the observational constraints from the luminosity 
function. 

Our model requires a transmission fraction (Tigm = 
0.5) higher than the favored values for a symmetric Lyman-a 
line (Tigm < 0.4). Due to the degeneracy between the escape 
fraction and the IGM transmission (the observed Lyman-a 
luminosity is reduced by a factor fa x Tigm) lower transmis- 
sion fractions Tigm can be compensated by higher escape 



fractions fa . Given that the dumpiness and amount of dust 
is already fixed by the LBG modelling, the freedom of chang- 
ing /„ is restricted. The most probable interpretation is that 
outflows have to be considered in order to provide a physi- 
cal picture that explains the values of Tigm ~ 0.5 we need. 
As mentioned in the last section, outflows will also have an 
impact on the fa escape fraction. Quantifying the effect of 
outflows is an issue we will addres in the future using the 
tools and simulations presented in Paper I and Paper II. 

At redshifts z > 6, before reionization is completed, 
the estimations of the IGM transmission fraction depend 
on the assumptions made to model both the reionization 
and t he Lyman-a line radiative transfer (|McQuinn et al] 
I2OO7I : iDaval et al] I2OIII : iDiikstra etal] l20l j ). This makes 
even more difficult to reach a conclusive prediction on the 
expected values for Tigm at a flxed redshift and will influence 
the interpretation of our results beyond z > 6, as will be 
discussed in Section 4.2. 



2.6 Comparison with the observational 
Luminosity Functions 

In Paper I we applied the above mentioned extinction model 
to match the rest-frame UV LPs derived from the simula- 
tion to the observed LPs in the redshift range 5 < z < 7. 
The flnal model imposes additional extinction in all the stel- 
lar populations younger than tc = 25 Myr, with parameters 
fj, — 0.01 for redshifts z ~ 5, 6 and /j, = 0.03 for redshift 
z ~ 7, and a = 1.5. At the faint end (Muv < —20) we 
have spotted a slight overabundance at all simulated red- 
shifts that could hint towards a low efficiency in the super- 
novae feedback that can regulate the star formation process 
in the shallow potential wells of the corresponding halos with 
masses ~ 2 x lO^°/i"^M0. 

In Paper II, with the dust contents already constrained 
from the UV LPs, we applied the model for the Lyman-a 
escape fraction to construct the LAE LPs. We note that 
the LAEs used to construct the LP had an equivalent width 
(EW) larger than a fiducial threshold of 20A, imposed to 
resemble a minimal threshold for observational detection. 
At redshifts 2 ~ 5 and 2 ~ 6 we found a good match at 
the bright-end within the Poissonian and cosmic variance 
uncertainties. The match with observations is also good if 
we consider IGM tr ansmission fraction s of 58 and 51 per 
cent as calcula ted bv lStark" et alJ (|201ll l using the results by 
iMeiksinI l|2006h . 

However, there is an excess at the faint end of the 
simulated LP with respect to the observations. The over- 
abundance corresponds to luminosities of Lbya ~ 2 x lO'''^ 
erg/s, where we found as well a mild over-abundance at the 
faint end of the UV LPs, if the star formation is lowered in 
some of these numerical galaxies as to match the observed 
abundance, both the intrinsic UV and Lyman-a emission 
will drop, making the galaxies to fall below the detectable 
range. Therefore, this change will not bear strong conse- 
quences for the fraction Xhya because of the simultaneous 
changes in LBGs and LAEs abundance. 

At redshift 2 ~ 7 the normalization of the LAE LP is 
still higher than observed at all luminosities by a factor of 
~ 0.4 in luminosity. By dimming each LAE in the simulation 
by a factor of /t ~ 0.4 we can repr oduce the LP abundances 
at 2 ~ 6.5 (jKashikawa et aLlbOOd ). 
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Figure 1. Fraction of galaxies showing strong Lyman-« emis- 
sion as a function of the absolute Myv magnitude in the mod- 
els (empty symbols) and observations (filled circles). The error 
bars in the model are calculated from the Poissonian variance on 
the number of strong emitting LAEs. The results of our model 
(squares) in the redshift range 5 < ^ < 6 shows a broad agree- 
ment with the observations reported by [stark et al 



3 RESULTS 

3.1 Lyman-a fraction at 5 < ^ < 7 

For each well resolved galaxy in the three simulation snap- 
shots at z = 5, 6, 7 we calculate the UV magnitudes cor- 
rected by dust absorption, Muv, the slopes of the spectra 
between 1200 A and 1600 A, (3, together with the intrinsic 
and observed Lyman-a luminosities. Then, for each galaxy 
we calculate the continuum flux red-wards of the Lyman-a 
line at 1240 A from the /3 slope values and the flux calculated 
from the Muv magnitudes. This way of calculat ing the con- 
tinuum was chosen to mimic the method used in lstark et al.l 
l|2010l ). The Lyman -a EW is then calculated as the ratio of 
the observed intensity in the Lyman-a line and the UV con- 
tinuum. 

We then bin the galaxies in Muv and calculate the frac- 
tion of LBGs with Lyman-a EW larger than 50 A. In Fig[T] 
we plot the fraction of galaxies with strong Lyman-a emis- 
sion (xhya) at redshifts z = 5 and z = 6 combined, and 
compared to the observational results at 4.5 < z < 6.0. As 
can be clearly seen in the plot, we are able to reproduce the 
observational trend without any flne-tuning of our simulated 
galaxies. Namely the values of Xhya for the bright LBGs with 
Muv < —20.5 and the increase in the fraction of LAEs for 
faint LBGs with Afuv > —20.5. The upturn in XLya around 
Muv = —21.0 is also reproduced, which is a consequence 
of reproducing the shape of the UV LP which shows the 
exponential drop in the abundance of LBGs beyond this 
magnitude, making the fraction Xhya at the brightest bin 
sensitive to rise at large values even for a low number of 
detected LBGs with strong Lyman-a emission. On aver- 
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Figure 2. Same as Figure [T] Filled symbols represent the ob- 
servational results by Schenker et al. (2011) for LBGs at z > 6.3 
with Lyman-o emission with EW> 25A. The empty symbols rep- 
resent our results from the simulation at z ~ 7. The theoretical 
data-points have been shifted by ±0.1 in redshift for clarity. The 
first scenario (circles) corresponds to the results of our model that 
do not include any evolution in the absorption model at z 7 
(/x = 1.0), while the second (triangles) adds a further dimming 
each of LAE by a factor of /t = 0.4. The factor was chosen to 
bring ou r results into agreemen t with the observed LAE LF at 
z ~ 6.5 l|Kashikawa et al.ll2006h . This scenario provides an im- 
proved match with the results reported by Schenker et al. (2011). 



age, a fraction x^ya — 0.47 ± 0.01 of the galaxies between 
—20.5 < Afuv < —18.5 show strong Lyman-a emission 
with EW> 50A, where only statistical (poissonian) errors 
have been considered in this estimate. 



3.2 Evolution at z > 6 

Recently, results of new line searches with deep Keck spec- 
troscopy of 19 LBGs in the redshift range 6.3 < z < 
8.8 selected with W FC3/IR data have been published 
llSchenker et al ]|201ll ). The spectroscopic exposures were de- 
signed to reach EW lower than 50 A. Under these condi- 
tions, only 2 LAEs are convincingly detected. These results 
suggest a strong evolution with respect to the observations 
at 5 < 2 < 6. Using complementary observational results 
from ultradeep optical spec troscopy with F0R S2 on the 
Very Large Telescope (VLT) (jFontana et al.ll2010 l). together 
with the results already obtained at z ~ 6 (jstark ct al] 
l201ll ). the authors derive a fraction of LAEs with EW> 25A: 
XLyc, = O.OSSto °f (-21.75 < Muv < -20.25) and XLyc = 
0.24±0.15 (-20.25 < A'/uv < 18.75). It is important to keep 
in mind that these fractions (estimated from the measure- 
ments at 6.3 < z < 8.8) are based on an extrapolation of 
the EW distributions at redshift z = 6. 

In Fig.[5]we present the results of our model for the x^ya 
fraction (EW> 25A) from the calculations with /t = 1.0 
(i.e. no evolution from the 2 = 6 extinction model) and 
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/t = 0.4. We recall (see the final paragraph of §2.6 ) that 
the factor /t is introduced to tune the spatial abundance 
of LAEs at 2 = 7 with respect to the LF inferred from 
observations 

Both models (/t ~ 1 and /t = 0.4) yield a simi- 
lar trend with Muv as the one found in 5 < z < 6. The 
fraction of LBGs with large EW increases for LBGs. How- 
ever, the model that matches the LAE LF at z ~ 7 with 
/t = 0.4 provides results for the xlvc fraction closer t he the 
results derived from observations bv ISchenker et all (|201ll ) 
with a;Lya = 0.03 ± 0.03 (-21.75 < Muv < -20.25) and 
XLyc = 0.37 ± 0.01 (-20.25 < Muv < 18.75). 



4 DISCUSSION 

In the previous sections we have presented a model for LBGs 
and LAEs that reproduce the observed trends for the x^ya 
fraction in two important aspects: 1) the rise of XLya for 
fainter Muv magnitude at a fixed redshift and 2) the overall 
evolution as a function of redshift between 5 < z < 7. In this 
Section we discuss the key elements in our model behind 
these results. 



4.1 Dependence of x^ya on Muv 

Assuming that the UV continuum and the Lyman-a come 
from star formation, one can express the intrinsic lumi- 
nosities as a function of the star formation rate: L^uv = 
1.4 X 10*° X SFR/(M0yr-l) erg and Ll yc = 

1.9 X 10"^ X SFR/Moyr-i erg s"^ (|Kennicutt|[l998l ). The 
exact conversion factors depend on the IMF and the metal- 
licity of the stellar population. We will take these values as 
constant, although this approximation is not fundamental 
in the argument that follows. 

Under these assumptions, the intrinsic equivalent width 
(EWi) is constant for all galaxies, EWi ~ 135A, making all 
bright LBGs detectable as strong LAEs regardless of Afuv if 
one applies the EW cut > 50A. The fact that only a fraction 
of bright LBG are also strong LAEs can be understood as 
a variation of the observed EW values from the intrinsic 
EWi. This is naturally expected in the presence of extinction 
and/or preferential scattering of Lyman-a photons out of the 
line of sight. 

We can assume that extinction reduces the intensity of 
the line by a factor fa and the UV continuum by a factor 
fc- Further dimming by the neutral IGM can also reduce the 
Lyman-a line intensity by a factor Tigm. In this case, the 
measured equivalent width is EW= Tigm x (/Q//c)xEWi. In 
our model we have taken Tigm to be constant for all galaxies 
at a given redshift. Under this assumption the properties of 
the EW distribution can be attributed to the scatter in the 
extinction factors fc and fa Q. 

Observationally, the factor fc can be inferred from the 
UV spectral slope, p. In the observed LBGs at z > 5 
l|Bouwens et al.ll2010l : iDunlop et a"l]|201ll ) and in our model 
as well, the evolution of the mean /3 slope at a given magni- 
tude Muv is not strong, despite the presence of considerable 



^ Although extreme values of the EW can be used to infer the 
presence of unusual stellar populations l lDiikstra &: WestrallioTol) 



scatter. On the other hand, the factor fa can vary at least 
by a factor of 10 depending on the properties of the galaxy 
(|Forero- Romero et al.l 120111 ). This suggests that, at a fixed 
redshift, the key physical factor in our model accounting 
for the trend of iLya with Afuv is the dependence of the 
Lyman-a escape fraction on galaxy luminosity. 

To test this assumption, we show in the left panel of Fig. 
[3]the results of approximating the Lyman-a escape fraction, 
/q , as a constant value. To construct this figure we keep the 
results for the UV continuum fixed and reduce the Lyman- 
a emission for all galaxies by a constant value fa ~ 0.2, 
a choice that also brings the LF for LAEs into agreement 
with observations at z ~ 5 and z ~ 6. This model fails to 
reproduce the observational trends for x^ya, in particular 
it seriously underestimates the fraction of bright LAEs for 
faint LBGs to be xi^ya = 0.10 ± 0.01. The Lyman-a escape 
fraction is too low for faint LBG, reducing too much the 
equivalent width and making most of these galaxies unde- 
tectable as LAEs. 

To further illustrate our point, we have also imple- 
mented a phenomenological model where the fa escape frac- 
tion is calculated from an effective optical depth to Lyman- 
a line photons, which is assumed to be proportional to the 
continuum optical depth. 



1 - cxp(-7Td) 



(6) 



where 7 > 1 and Td is the continuum o ptical depth at a wave- 
length of 1216A. [Shimizu et al.l (|201ll ) have recently applied 
this model to reproduce the properties of LAEs at z = 3.1. 
In Fig. [3] we show the results for 7 = 50, this value was 
chosen to give a close match to the LAE LF at 2 ~ 5 and 
z ~ 6. This model includes a dependence of the fa escape 
fraction on galaxy gaseous mass (and hence luminosity) , via 
the dependence on the gas optical depth that is larger for 
luminous systems. 

This approximation manages to reproduce similar re- 
sults for xi^ya as our radiative transfer motivated model, a 
result that is not surprising given the fact that the functional 
form of Eq.® is close to the expected from radiative trans- 
fer effects of the Lyman-a line in the presence of neutra l 
Hydrogen (|Hansen fc Ohll2006l : Forero-Romero et aLlbOllI ). 

The result from this phenomenological model adds sup- 
port to our claim that the key factor in our model behind 
the shape of the a;Lya-Muv plot is the mass dependence of 
the Lyman-a escape fraction that makes luminous galaxies 
have, on average, small values < 0.2 while fainter galaxies 
have large fa > 0.2. 



4.2 The evolution of XLya beyond z > 6 

Different authors report a change in the properties of LAE 
emitting galaxies b eyond redshift z > 6, both in the LF 
(|Haves et al. (I2OII) and refe rences therein) and the xi^ya 
fraction ( Schenker et al.|[2011^ . 

In our model, the d rop in the LAE abundance at z ~ 6.5 
2OO6I I and the evolution of Xi^ya at z > 6.3 



Kashikawa et al.l. 

Schenker et al.ll201lh re explained by the evolution in the 
IGM transmission fraction. Such evolution in the transfer 
of the Lyman-a line is required in our model to match the 
observational constraints. The dimming of all LAEs at z ~ 7 
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Figure 3. Same as Fig. [T] for two different models of the Lyman-a escape fraction. In both cases we keep the results for the Muv 
magnitudes fixed and vary the Lyman-a escape fraction. The upper panel show s the results of a co nstant escape fraction fa = 0.2, which 
clearly present a trend in contradiction with the observational constraints from lStark et ahl ||2010| ). The lower panel shows the results of 
the model based on radiative transfer results (squares) compared against a phenomenological model where the Lyman-o escape fraction 
is calculated by Eq|6] (circles). The phenomenological approach shows a trend similar to the one obtained the radiative transfer model. 



by a constant factor /t = 0.4 makes the results of our model 
to better fit the LAE LF at z ~ 6.5 and the xi^ya fraction 
at z > 6.3. 

If confirmed, the observational evidence for a decrease 
in the observed XLya fraction at redshifts z > 6 can be 
attributed to three reasons: 

1) an increase in the extinction in galaxies ISM 

2) evolution in the Hydrogen neutral fraction in the IGM 

3) an increase in the escape fraction of ionizing photons. 



Given the observational and t heoretical results for 
UV continuum slopes for z > 6 (iBouwens et al.1 |2010| : 
iDunlop et alll201ll : iForero-Romero et ahlbOloF it is rather 
implausible that the extinction has increased. Therefore, 
only the explanations 2) and 3) seem more plausible. It is 
beyond scope of this paper to distinguish between these two 
scenarios. 

However, if we add up the fraction /t to the full trans- 
mission effect through a neutral IGM we have 7igm(2 ~ 
7.0) = /t X 0.5 — 0.2, where 0.5 corresponds to the fiducial 
value for the IGM transmission used in our model. The in- 
terpretation of this fraction Tigm = 0.2 in terms of a global 
neutral fraction Xhi is considerably more difficult. Different 
models can give results in the range Xhi ~ 0.4—0.9 depend- 
ing on the treatment of the line through the IGM, making 
it difficult to rule out large values of the neutral HI fraction 
from the evolution in XLya, specially when gas bulk veloci- 
ties f r om winds are includ ed in the analysis (jMcQuinn et al.l 
I2OO7I : iDiikstra et aLllioill ). 



5 CONCLUSIONS 

Observations report that not all bright Lyman Break Galax- 
ies in the redshift range 5 < z < 7 are seen as strong Ly- 



man Alpha Emitters ("St ark et al.]|2010l . l201ll : lFontana et al.l 
I2OIOI : iSchcnkcr ct al. 201ll ). The fraction of LBGs that can 
be detected with strong Lyman-a emission, xi^ya tests the 
equivalent width distribution, providing additional informa- 
tion that is not included in the luminosity function of LAEs 
and LBGs. In this paper we present a theoretical model of 
high redshift LBGs and LAEs that reproduce the observed 
trends for xi^yc, in the redshift range 5 < z < 7. We find that 
the a decreasing Lyman-a escape fraction with increasing 
galaxy luminosity is a key element in our model to explain 
the observations. 

In the redshift range 5 < z < 6 observations show a 
clear evolution of the x^ya with absolute magnitude Afuv- 
From the results of our model, we suggest that a key ingredi- 
ent to explain this trend is the decrease of the Lyman-a es- 
cape fraction with increasing galaxy luminosity. We test this 
hypothesis by fixing the UV continuum results and assuming 
a constant value for . This produces an increasing fraction 
of LBGs with strong Lyman-a emission for brighter Afuv, 
the exact opposite trend of what is reported by observations. 
In contrast, a model with a decreasing Lyman-a escape frac- 
tion with increasing galaxy mass can simultaneously repro- 
duce the observed trends for Xhya and the LAE-LF. The ad- 
ditional influence of an IGM transmiss ion coefficient, Tig m, 
dependent on large scale environment ijZheng et al.ll201ll ) is 
a point that deserves further investigation. In particular, it 
would be useful to know to what extent this kind of envi- 
ronment dependent transmission explains the trend of x^ya 
with absolute magnitude Afuv- 

In order to reproduce the xi^ya evolution at z > 6.3 
our model requires a decrease of the transmission fraction 
of Lyman-a photons by a factor of fr ~ 0.4, also improv- 
ing the agreement with the observed LAE LF at z ~ 6.5, 
although confirmation of such evolution from increased sen- 
sitivity observations and larger samples is still needed. 
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The fraction xi^ya provides an additional statistical tool 
to describe galaxy populations evolving during the reioniza- 
tion epoch. More sensitive instruments and larger galaxy 
samples will be necessary to quantify the evolution of the 
XLyct fraction. In the same spirit, this statistic can be used 
as a benchmark for a theoretical model aiming at describing 
the evolving galaxy populations at these redshifts. 
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